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Abstract 

The influence of a Zeeman magnetic field on the superconducting char- 
acteristics of the attractive Hubbard model was investigated. The ground 
state and temperature phase diagrams were obtained for a fixed number 
of particles. Two critical magnetic fields were found for the first order 
phase transition from the superconducting to the normal state for n 7^ 1. 
For some range of parameters a reentrant transition was found and gap- 
less superconductivity was obtained. The significance of the Hartree term 
was also analyzed. 

1 Introduction 

The impact of the Zeeman coupling between the spins of the electrons and an 
applied magnetic field on superconductivity has been analyzed for many years 
[l]-[5]. One of the well-known results concerning this influence is the discovery 
of the so-called Clogston limit. In the weak coupling regime (s-wave pairing), at 
T = 0, the superconductivity is destroyed through the paramagnetic effect and 
there is a first-order phase transition to the normal state at a universal value of 
the critical magnetic field he — Ao/\/2 « 0.707Ao [1], where Aq is the gap at 
T = 0, ft. = 0. Further investigations have revealed that in the presence of mag- 
netic field the formation of Cooper pairs across the spin-split Fermi surface with 
non-zero total momentum (fc |, —k + q i) is possible, leading to the so-called 
Fulde and Ferrell [3], and Larkin and Ovchinnikov [4] (FFLO) state, which is 
favored (against the normal state) up to hf^^^ = 0.754Ao. Because of the very 
strong destroying influence of the orbital effect on the superconductivity, it is 
hard to observe the above-mentioned effects in ordinary superconductors. How- 
ever, the development of experimental techniques in cold atomic Fermi gases 
allows investigation of the spin polarization on superfluidity [6]- [11] in a sim- 
ple way. Recent experiments [11] have indicated that in the density profiles of 



1 



trapped Fermi mixtures with population imbalance there is an unpolarized su- 
perfluid core in the center of the trap and a polarized normal state surrounding 
this core (phase separation). 

In this paper we analyze the influence of the Zeeman term on the super- 
fluid characteristics of a lattice fermion (the spin-polarized attractive {U < 0) 
Hubbard) model [12]: 

H = Yl^^ij ~ ^^y)4<rCja + UY.nitnn - hJ2{nit - "ii), (1) 

ija i i 

where: a =t, i, nia = cl^Cia, t^j - hopping integral, U - on-site interaction, fj, - 
chemical potential. The Zeeman term can be created by an external magnetic 
field (in (g^Bf^) units) or by a spin population imbalance in the context of cold 
atomic Fermi gases. 

We consider the case of pairing only with q = 0. The gap parameter is 
defined by: A = -^Y.,{c^lC,^) = --^ Efe(c_fe|%)- Applying the broken 
symmetry Hartree approximation, we obtain the grand canonical potential f2 
[13] and the free energy F: 

I = \un{2 - n) - m(1 - n) + \uM^ - ^ (2) 

fe ^ ' 

where: N mmiber of the lattice sites, f3 — l/ksT, Af = rif — rij - spin magne- 
tization (polarization), n^, = jj X]fe('^^ '^^o-) spin- up and spin-down electron 
density, n — n]+ni - electron concentration, Ej:^ = (— f^+t^)8j:-|-^^-|-/i-|-a;g, 
E^^^ = (-tT + t^)Q. _ m _ /, + ^- ^ v'((~^T_ti)eg-A)2 + |A|2 are 

quasiparticle energies, p = /j. — Here 9^ = J2f=i cos{kiai) (d = 3 for three- 
dimensional lattice), a; = 1 in further considerations. Using (2), one can easily 
get the equations for the gap, particle number (determining jj) and magnetiza- 
tion [13]: f£ = 0, ^ = 0, M = — -i-^, respectively. In the above equations, 
we have taken into account the spin dependent Hartree term. This method can 
be called a BCS-Stoner approach. In the following, we set = = t and use 
t as the unit. 

2 Results 

We have performed an analysis of the influence of magnetic field on superfiuidity, 
focusing on the ground state and the temperature phase diagrams for a fixed 
electron concentration and indicating differences with respect to the case of a 
fixed chemical potential [13]. Only selected results will be presented below, 
while an extended paper will be published elsewhere [14]. 
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Figure 1: The critical magnetic field vs. the chemical potential (a) and the 
electron concentration (b) for the first order SC-NO transition, at T = 0; NO-I 
- partially polarized (P < 1) and NO-II - fully polarized {P = 1) normal state, 
e - empty state, PS - phase separation. In (a) solid lines denote SC-NO first 
order transition, dashed lines border NO-I and NO-II states, borders between 
NO-II and empty states are shown by dash-dotted line. Charge density wave 
(CDW) state, being degenerated with SC for h=0, n=l, is not shown. 



Fig. 1 shows the ground state phase diagrams for a simple cubic lattice, 
i.e. critical magnetic fields for the first order superconducting-normal (SC-NO) 
phase transition vs. the chemical potential (a) and the particle concentration 
(b), for three values of the on-site attraction, ranging from a weak to interme- 
diate coupling. If the number of particles is fixed and n ^ 1, one obtains two 
critical magnetic fields in the phase diagrams [13]. The two critical fields (Fig. 
lb) define the phase separation (PS) region between the superconducting phase 
with the particle density and the normal state with the density of particles 
n„, as opposed to the case of a fixed chemical potential (Fig. la). Moreover, 
the increasing on-site attraction widens the range of occurrence of PS. One can 
also distinguish the partially polarized (P = (n^ — n|)/(n| + n^) < 1) nor- 
mal state (NO-I) and the fully polarized (P = 1) normal state (NO-II) in the 
phase diagrams. The same regions appear for n > 1 because of the particle-hole 
symmetry. The first order transition lines at P = have been determined nu- 
merically from the condition fl^^'^ = ^n^^ (where and fl^^^ denote the 
grand canonical potential of the normal (A = 0, P ^ 0) and the superconduct- 
ing (A ^ 0, P = 0) state, respectively.). 

Fig. 2 shows the magnetic field-temperature {h — T) phase diagram (3D case) 
for the fixed n — 0.75 and U = —4. The lower set of curves denotes the diagram 
without the Hartree term. With increasing magnetic field, the character of the 
transition between SC and NO state changes from the second order (solid line) 
to the first order, which starts from the tricritical point (TCP). We can also 
observe the region of the phase separation between SC (A 7^ 0, P 7^ 0) and 
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NO (A = 0, P ^ 0) state in contrast to the (T — h) phase diagrams obtained 
for a fixed chemical potential [13]. The presence of the Hartree term raises the 
values of the critical magnetic field at T = for the first order transition, which 
exceed the Clogston limit. In the phase diagram without the Hartree term 
(conventional BCS like) the PS region is very narrow for the chosen values of 
the parameters. For small \U\ and n or for n = 1, the Clogston limit is achieved 
in the phase diagrams without the Hartree term. For sufficiently high fields a 
reentrant transition takes place in the phase diagram with the Hartree term. 
Hence, the increase in temperature can induce superconductivity. We can infer 
that the Hartree term can be of importance in the context of the spin-polarized 
superconductivity and the neglect of this term may lead not only to quantitative 
but also to qualitative changes in the phase diagrams. We have also found the 
gapless superconducting (GSC) region (Fig. 2) i.e. spin-polarized SC, which 
has a gapless spectrum for the majority spin species for some magnetic fields 
{h > A) and temperatures. 




Figure 2: Magnetic field-temperature phase diagrams for fixed n — 0.75, U = 
—4, for a simple cubic lattice. The lower set of curves denotes the diagram 
without the Hartree term. The arrow indicates TCP. The thick dotted lines 
denote the first order phase transition. These lines bound PS region. GSC - 
gapless superconductivity. The thin dotted line is merely an extension of the 
2nd Qj-jgj. transition line below TCP (metastability line). 
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3 Conclusions 



We have considered the influence of a Zeeman magnetic fleld on superfluid prop- 
erties of the attractive Hubbard model. The analysis has been restricted to the 
case of s-wave pairing with q = 0. For a fixed number of particles and n ^ 1, 
one obtains two critical magnetic fields in the phase diagrams, which limit the 
PS. In the PS region not only the particle densities but also the polarizations of 
the coexisting SC and NO states are different. We also observe a reentrant tran- 
sition for sufficiently high magnetic fields in the h — T phase diagram with the 
Hartree term. Moreover, the presence of the Hartree term causes an increase of 
he of the Clogston limit. At T > the spin-polarized SC state and the gapless 
superconducting region have been found in h — T phase diagrams. 

Acknowledgements 

We would like to thank Prof. S. Robaszkiewicz for discussions. 

References 

[1] A. M. Clogston, Phys. Rev. Lett. 9, 266 (1962). 

[2] G. Sarma, J. Phys. Chem. Solids 24, 1029 (1963). 

[3] P. Fulde, R. A. Ferrell, Phys. Rev. 135, A550 (1964). 

[4] A. I. Larkin, Y. N. Ovchinnikov, Zh. Eksp. Teor. Fiz. 47, 1136 (1964). 

[5] L. W. Gruenberg, L. Gunther, Phys. Rev. Lett. 16, 996 (1966). 

[6] D. E. Sheehy, L. Radzihovsky, Phys. Rev. Lett 96, 060401 (2006). 

[7] K. B. Gubbels, M. W. J. Romans, H. T. C. Stoof, Phys. Rev. Lett 97, 
210402 (2006). 

[8] Ch. Chien, Q. Chen, Y. He, K. Levin, Phys. Rev. Lett. 98, 110404 (2007). 

[9] M. M. Parish, F. M. Marchetti, A. Lamacraft, B. D. Simons, Phys. Rev. 
Lett 98, 160402 (2007). 

[10] T. K. Koponcn, T. Paananen, J. P. Martikainen, P. Torma, Phys. Rev. Lett 
99, 120403 (2007). 

[11] Y. Shin, M. W. Zwierlein, C. H. Schunck, A. Schirotzek and W. Ketterle, 
Phys. Rev. Lett. 97, 030401 (2006); Y. Shin, C. H. Schunck, A. Schirotzek, 
W. Ketterle, Nature 451, 689 (2008). 

[12] R. Micnas, J. Ranninger, S. Robaszkiewicz, Rev. Mod. Phys. 62, 113 (1990). 



5 



[13] A. Kujawa, R. Micnas, Acta Phys. Pol. A 114, 43 (2008). 
[14] A. Kujawa, R. Micnas, in preparation. 



6 



